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Abstract:  Stepped  aperiodic  waveguide  structures  are  designed  for  mode  conversion  by 
numerical  optimization  coupled  with  a  mode  matching  solver.  Several  configurations  and 
design  goals  are  evaluated  for  performance  and  dimensional  tolerances.  Results  are  validated 
against  a  FDTD  code. 

Introduction 

Previous  studies  by  Haq,  Webb  and  others  have  demonstrated  the  design  of  aperiodic 
waveguide  structures  to  act  as  Liters  and  mode  converters  [1,  2],  These  aperiodic  structures 
have  been  shown  to  yield  high  efficiency  mode  conversion  or  filtering  in  lengths  considerably 
shorter  than  structures  using  gradual  transitions  and  periodic  perturbations.  The  design 
method  developed  by  Haq  and  others  has  used  mode-matching  models  [3]  for  the  irregular, 
stepped  waveguides  coupled  with  computer  optimization  to  achieve  the  design  goal  using  a 
Matlab  optimization  routine. 

Similar  designs  are  described  here,  using  a  mode  matching  code  written  in  Fortran  and 
with  optimization  accomplished  with  the  downhill  simplex  method  with  simulated  annealing 
using  an  algorithm  from  the  book  Numerical  Recipes  in  Fortran  [4],  Where  Haq  et  al.  looked 
mainly  for  waveguide  shapes  with  relatively  wide  cavities,  we  have  sought  lower  profile 
designs.  It  is  found  that  lower  profiles  can  meet  the  design  goals  and  result  in  a  structure 
with  lower  Q.  In  any  case,  there  appear  to  be  very  many  possible  configurations  for  a  given 
mode  conversion  goal,  to  the  point  that  it  is  unlikely  to  find  the  same  design  twice.  Tolerance 
analysis  was  carried  out  for  the  designs  to  show  edge  sensitivity  and  Monte  Carlo  degradation 
rate.  The  mode  matching  code  and  mode  conversion  designs  were  validated  by  comparison 
with  FDTD  solutions  for  the  discontinuous  waveguides. 

Optimization  for  design  of  waveguide  structures 

In  the  mode  matching  solution  the  fields  in  two  waveguides  forming  a  junction  are 
expanded  in  a  sum  of  modes  and  the  continuity  of  E  and  H  fields  is  enforced  as  integrals 
weighted  with  the  mode  functions.  The  procedure  described  in  [3]  yields  the  scattering 
matrix  for  a  junction.  The  scattering  matrix  for  a  sequence  of  junctions,  taking  account 
of  all  interactions,  can  then  be  obtained  by  cascading  individual  scattering  matrices.  We 
found  that  rough  results  could  be  obtained  with  modes  TEio  and  TE30  in  the  input  guide, 
where  only  TE10  was  above  cutoff,  and  the  number  of  modes  in  other  sections  scaled  in 
proportion  to  width.  Higher  accuracy  was  obtained  with  modes  up  to  TE90  in  the  input 
and  proportional  scaling.  The  mode  matching  solution  took  about  0.03  to  0.1  seconds  for  a 
single  frequency  with  20  waveguide  junctions  on  a  1.7  GHz  Pentium  computer  using  modes 
up  to  TE90  in  the  input.  These  wave  guide  structures  were  also  modeled  with  a  2D  FDTD 
code  with  a  plane  of  symmetry.  The  FDTD  solution  could  take  20  minutes  to  several  hours 
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to  yield  all  frequencies  for  a  given  waveguide  configuration,  but  could  offer  options  for  more 
generality  in  the  model  if  needed.  Both  mode  matching  and  FDTD  codes  were  conpled  with 
optimization,  but  the  mode  matching  is  of  course  much  nicer  to  use. 

To  synthesize  waveguide  mode  converters  the  solver  was  coupled  with  an  optimization 
code  that  used  the  downhill  simplex  method  with  simulated  annealing  [4].  For  an  N  pa¬ 
rameter  space,  a  simplex  is  a  figure  of  N  +  1  vertices  and  connecting  edges,  in  this  case 
N  +  1  different  configurations  of  the  waveguide  sections.  The  objective  function  is  evaluated 
at  each  vertex.  Then  the  algorithm  may  reflect  the  highest  (worst)  vertex  relative  to  the 
others,  and  may  also  expand  or  contract  the  simplex.  The  initial  simplex  was  generated  by 
entering  an  initial  waveguide  profile  and  generating  N  additional  profiles  with  random  per¬ 
turbations.  Simulated  annealing  adds  a  randomness  to  simulate  a  metal  cooling  to  its  lowest 
energy  state.  It  is  stated  in  [4]  that  an  annealing  schedule  must  be  chosen  for  each  particular 
problem.  We  usually  just  optimized  with  zero  temperature  from  the  initial  random  simplex, 
and  if  the  result  looked  promising,  then  raised  the  temperature  and  lowered  it  to  see  if  a 
better  solution  could  be  found. 

Two  mode  converter  designs  are  shown  in  Fig.  1  through  4  for  a  goal  to  optimize  TEio 
output  at  9.6,  9.7  and  9.8  GHz  and  TE30  output  at  10  GHz  for  TE10  input.  The  cost 
function  was  the  product  of  the  output  powers  for  the  designated  modes  and  frequencies. 
The  waveguide  had  20  sections  of  1  cm  length  that  were  only  adjusted  in  width.  Both  designs 
meet  the  goals  fairly  well,  although  the  second  has  a  lower  tolerance.  The  result  of  a  FDTD 
solution  for  the  final  design  is  included  in  Fig.  1,  and  is  in  good  agreement  with  the  mode 
matching  result.  A  design  to  pass  TE10  at  9.7  and  10.3  GHz  and  convert  to  TE30  at  10  GHz 
is  shown  in  Fig.  5  and  6.  In  this  case  the  waveguide  had  10  sections  that  were  adjusted  in 
both  width  and  length.  For  this  goal  we  were  unable  to  find  a  good  low  profile  design  with 
20  sections  adjusted  in  width.  We  did  End  a  good  solution  with  21  cm  maximum  width,  but 
it  had  a  very  low  tolerance.  A  good  low  profile  design  was  also  found  with  40  sections  of  0.5 
cm  length  adjusted  in  width. 

Conclusion 

The  mode  matching  solution  for  discontinuous  waveguides,  coupled  with  the  downhill 
simplex  method  of  optimization  and  simulated  annealing  provides  an  efficient  means  of  de¬ 
signing  waveguide  structures  as  mode  converters  and  filters.  In  many  cases  good  results  were 
obtained  in  very  reasonable  computer  time  by  using  a  minimal  number  of  modes.  When  a 
promising  result  was  obtained  the  number  of  modes  could  be  increased  to  ensure  accuracy. 
A  similar  strategy  was  described  in  [2].  There  may  be  many  possible  solutions  for  a  given 
goal  and  many  not  so  good  relative  minima,  but  the  method  of  Simulated  Annealing  was 
helpful  in  finding  a  good  design. 
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Fig.  1.  Waveguide  mode  converter  optimized  for  TEi0  output  at  9.6,  9.7  and  9.8  GHz  and  for  TE30  output 
at  10  GHz  for  TE10  in  the  input  guide.  The  optimization  was  started  with  a  flat  initial  profile  with  width  of 
8  cm  and  20  adjustable  sections  with  1  cm  lengths.  The  input  guide  on  the  left  has  width  2.29  cm,  and  the 
output  guide  has  width  5.97  cm.  The  “2a;”  scale  shows  the  full  width  of  the  waveguide. 
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Fig.  2.  Dimension  tolerance  plots  for  the  waveguide  in  Fig.  1,  including  Monte  Carlo  analysis  for  er¬ 
ror  range  ±5  and  edge  sensitivity  A  for  one  percent  degradation  in  the  cost  function  C  where  C  = 
[Hteio(9-6)Pteio(9-7)Pteio(9.8)Pte3o(10)]1^4  and  PrEno(/)  indicates  the  normalized  power  in  mode  TE„0 
at  f  GHz. 
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Fig.  3.  Waveguide  mode  converter  optimized  for  TEi0  output  at  9.6,  9.7  and  9.8  GHz  and  for  TE30  output 
at  10  GHz  for  TE10  in  the  input  guide.  The  optimization  was  started  with  a  tapered  initial  profile. 
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Fig.  4.  Dimension  tolerance  plots  for  the  waveguide  in  Fig.  3,  including  Monte  Carlo  analysis  for  error 
range  ±<5  and  edge  sensitivity  A  for  one  percent  degradation  in  the  cost  function  C.  TEno(/)  indicates  the 
output  power  in  mode  TE„o  at  frequency  /  GHz. 
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Fig.  5.  Waveguide  mode  converter  optimized  for  TE10  output  at  9.7  and  10.3  GHz  and  TE30  output  at  10 
GHz  for  TE10  in  the  input  guide.  The  optimization  was  started  with  a  flat  initial  profile  with  width  of  8  cm, 
and  10  sections  with  both  width  and  length  varied. 
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Fig.  6  .  Dimension  tolerance  plots  for  the  waveguide  in  Fig.  5,  including  Monte  Carlo  analysis  for  error 
range  ±(5  and  edge  sensitivity  A  for  one  percent  degradation  in  the  cost  function  C . 


